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' We explore analogy between two-component quantum alkali gases and spin-polarized helium 

systems. Recent experiments in trapped gases are put into the frame of the existing theory for 
Castaing instability in transverse channel and formation of homogeneous precessing domains in 
spin-polarized systems. Analogous effects have already been observed in spin-polarized 3 He and 
, 3 He — 4 He mixtures systems. The threshold effect of the confining potential on the instability 

q_) ■ is analyzed. New experimental possibilities for observation of transverse instability in a trap are 

["T i | discussed. 

*n : 

i 1! I. INTRODUCTION 

Remarkable observation of spontaneous spatial separation of states in recent experiments on trapped 87 Kb Q was 
immediately followed by three, almost simultaneous theoretical explanations Though these explanations differ 

in details, all three exploit the analogy between a two-component s7 Rb gas and a spin-polarized spin-1/2 gas. The 
analogy is based on the similarity of strong molecular fields in all two-component quantum gases in the Boltzmann 
temperature range. In spin-polarized gases, this molecular field can cause transverse spin currents Q which, in turn, 
can result in spatial separation of spin-up and spin-down components. An analog of this effect was, according to Refs. 
[§-fI| , observed in experiment jl| with 87 Rb. 

The aim of this paper is to point out that the analogy between multi-component and spin-polarized gases is even 
richer. In mid-90s there were several interesting experimental and theoretical publications on spin-polarized systems 
that are similar in spirit to Refs. The key is the so-called Castaing instability of systems with inhomogeneous 

polarization with respect to transverse fluctuations and the formation of precessing magnetic domains. These effects 
, have been studied, both experimentally and theoretically, in the context of spin-polarized helium systems. We will 
■ analyze the analogs of these effects for alkali gases and make suggestions for new types of experiments including the 
observation of the instabilities. So far, there have been no such instability experiments for alkali gases. Though the 
£S| ■ two-component alkali gases are to a large extent similar to spin-polarized helium, the characteristic parameters are 
different and the issue should be revisited in application, for example, to 87 Rb. The main difference is the small size 
, and non-uniformity of confined clouds of alkali gases. 
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II. MAIN EQUATIONS 



The analogy between spin-polarized quantum gases and multi-component gases (gases of particles with discrete 
internal degrees of freedom) has been discussed in detail in Ref. ||. For example, the densities of components of a 
two-component gas (e.g., 87 Rb) play the role of the spin-up and spin-down components of the spin-polarized gas, while 
the (Rabi) transitions between the states are similar to the transitions between the spin components and lead to the 
formation of the transverse (to the quantization axis) magnetization and correspond to the off-diagonal components 
of the spin density matrix. 

First, let us translate the equations of motion of Ref. |J onto the language of trapped two-component quantum 
Boltzmann gas in the s-wave approximation (in experimentally interesting case, the gas temperature is so low that 
the particle wavelength is much larger than the scattering length a). For a two-component gas, the density n is a 
2x2 matrix. The diagonal components ni,2 of this matrix represent densities of pure components and obey the 
conservation law 

9 t m )2 (r,t) + V i Jp ) (r,t) = 0, (1) 

( 1 2) 

where J\ ' is the corresponding current. The off-diagonal (mixed) matrix element riyi describes transitions between 
the components and obeys the equation 

d t n 12 (r,t) - i5Un 12 {v,t) + V< J. ( (12) (r, t) = Iff, (2) 
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where 

SU = E7i(r) - U 2 (r) + ( 5ll - ff22 )n(r), (3) 

Ui t 2 is the trapping potential for the component (1,2), and the interaction constants git in 87 Kb are expressed via 
the corresponding scattering lengths as 

gik = 47ra ife /m, Ou = o + a', a 2 2 = a - a', ai2 = a. (4) 

Since the asymmetry of scattering is weak, a' <C a, the collision integral /jjf^, which is proportional to (a'/a) 2 , 
is negligible, 7f| H ~ @. It is convenient to use, instead of matrix equations, the scalar and vector pseudo-spin 
components of the densities and currents as 

n = n\ + n 2 ,M =- (n 12 + n 2 i, in 2 i - in 12 ,nx - n 2 ) (5) 

and the same for the currents Ji and J^, respectively. In these pseudo-spin notations, the "magnetic field" or "Larmor 
frequency" has the form 

n L = (o,o,su). (6) 

In s7 Rb this field is low, £Il/T <C 1. Here, in contrast to experiments with true spins, it is almost impossible to 
change the direction of the magnetic field without exciting the Rabi transitions by rf pulses. 

In the hydrodynamic limit of small gradients, equations of motion ([!]) -(3) reduce to the well-known Leggett equation 
for the magnetic moment M(£, r) of the spin-polarized gas (see, e.g., ||3 and references therein): 

3 t M + Jl t xM + flA = 0, (7) 

d t J k + Sl L xJ k + 2gM xJ k + (T/m) j fc - -J fe /r, (8) 

j k = V fe M + MV k U/T + (n/4T) V fc O L (9) 

where U = (1/2) (U\ + U 2 + 3gn), g = Ana/m, and r is the relaxation time (in this case, the difference between 
transverse and longitudinal relaxation times is vanishingly small). In equilibrium, when the currents are absent, 
jfc = 0, and the "polarization" M follows the field and is directed along it: 

j fc 0) = V fe Mo + [MoVfef/ + n VfcO L /4] /T = 0. (10) 

In Eqs. (0)-@, the Larmor frequency can contain weak coordinate dependence, VLl = f^o + (r). Under usual 
assumptions, the terms dtJk x Jfc — 0, and Eq.(||) yields the following expression for the spin current: 

T ' ^T72 [jfc+^Jfc x Mi +/i 2 Mi (Mi .j fc )] , (11) 



1 + ii?M( 
IjM 1 = T(2gM+5n L ), 

where D = rT/m is the spin diffusion coefficient. 

If the magnetic field and trapping potential are uniform, Sflr, = 0, V k U = 0, Eq. ( |ll|) is exactly the same as the 
standard Leggett equation for the spin current. 



III. FORMATION OF DOMAINS 



Often, the magnetization in polarized helium can be non-uniform and noticeably different from its equilibrium 
field-driven value. The non-equilibrium polarization can be created by several techniques such as by flipping the 
magnetization by 180° by a proper rf pulse in a part of a cell, by optical pumping, by connecting two cells with 
different magnetizations, etc. This initially non-uniform and non-equilibrium polarization can rapidly evolve, by 
means of non-dissipative coherent spin currents, into a quasi-equilibrium state, which will then slowly relax towards 
the real equilibrium. Such a quasi-equilibrium state has been thoroughly studied theoretically and experimentally by 
Fomin and Dmitriev Q . They demonstrated that a large initial non-equilibrium polarization M (r) can lead to the 
coherent formation of longitudinal magnetization domains with spatially separated up and down spins. The effect was 
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observed in low temperature 3 He — 4 He mixtures and normal liquid 3 _ffe. The existence of such domains requires 
that the width of the domain wall A, 

A ' ~ M M|VO L |' (12) 

should be smaller than the system size, A <C L. Diffusive dissipative processes lead eventually to the homogenization 
of the system and the disappearance of domains. Numerical analysis the domain formation can be found also in Ref. 
fLOf . Note, that the field gradient is crucial for maintaining the quasi-stationary domain state; without this 

gradient, a domain- like distribution of magnetization becomes unstable with respect to transverse perturbations (see 
the next section). 

Though the main equations for spin-polarized helium and two-component alkali gases are similar, there are some 
differences in domain formation between helium and JILA experiments. In JILA experiment pi, the initial condition 
corresponds to the large transverse magnetization M tr ~ n/2 created by a 90° r/pulse. In experiments and calculations 
of Refs. 0, the initial state is longitudinal. However, since the formation of domains in Refs. || goes through 
large transverse currents, the quasi-equilibrium domain state is the same. A more important difference is that the 
experiments on alkali gases are performed only in small-size traps. Taking into account the trap profile U (x) ~ 
mu 2 x 2 /2 and that /i = 2gr, D = Tr/m, 5Hl = mv 2 x 2 /2, we get that in the trap of size x ~ L = y/T/mujQ, the 
width of the domain wall is 
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AnanL 2 v 2 
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(13) 



The typical values of the parameters, r ss 10 2 s, n w 4 • 10 13 cm 3 , a = 5 • 10 7 cm, m — 1.5 • 10 22 g, ujq ~ 7 Hz, 
combined with an estimate v < 0.1 Hz, lead to 

X/L > 0.34. (14) 

Therefore, it is conceivable that the separation of components in the JILA experiment is similar to the formation of 
three Fomin-Dmitriev domains. This picture of domain formation is consistent with the results of Refs. |J. The 
temporal symmetry of formation and disappearance of the domains is explained by the fact that in JILA experiments 
/j,M is not very large and the coherent processes, which are responsible for the formation of the domains, have 
approximately the same time constants as diffusion, which is responsible for the homogenization. Note, that because 
of the peculiar initial condition - large transverse polarization - the formation of domains does not require the instability 
in transverse channel as it is often necessary when starting from a non-equilibrium longitudinal configuration (see 
below) . 

There is another, earlier unexplained experimental fact. The increase of the field gradient v resulted in threshold 
stratification at some value v > 0.18 Hz Q|. The above estimates show that at v ~ 0.2 Hz the number of domains 
can jump to five and give the impression of such stratification. Since the diffusion equilibration accelerates rapidly 
with an increase in the number of domains, this can also produce an impression of instability in the experiment. 

IV. TRANSVERSE INSTABILITY: APPLICATION TO THE JILA EXPERIMENT 

It is worth noting that the dynamics of the domain formation in the experiment can often start as transverse 
instability in accordance with Castaing's mechanism |llj. Let us, first, briefly describe the Castaing instability in 
spin-polarized helium systems and then look at the implications for the JILA-type experiments. 

Suppose the system is homogeneously polarized with uniform magnetization Mo. If one is interested in small 
fluctuations <5M in transverse channel that are perpendicular to the equilibrium magnetization Mo, the linearized 
Eqs. (0), ([ll]) describe the spectrum of the circularly polarized Silin spin waves, 

Dk 2 

LU = Q ° + 1+/i 2 M 2 (* - ^0) ■ (15) 

When the magnetic moment is driven from the equilibrium and acquires a gradient VMo without change in direction, 
it is often assumed that this inhomogeneity will slowly disappear as a result of some relaxation process. This is not 
always the case. Often, the magnetic moment develops spontaneously a noticeable transverse component, which, in 
turn, may cause large longitudinal spin currents. Final relaxation towards the equilibrium is much slower. Castaing 
jllj showed that a purely longitudinal state with a non-equilibrium gradient of magnetization becomes unstable with 



3 



respect to transverse fluctuations (5M. The instability occurs with respect to longwave transverse fluctuations with 
the wave vectors k 2 < 2fik ■ VMq, Mathematically this means that the imaginary component of the spectrum (|l5| ) 



acquires, according to Eq. (11), the additional term 



The imaginary part of the spectrum of transverse fluctuations near the distribution M becomes negative when 

1 + ^M Z 



(17) 



The change in sign of Im u> signals the instability. Since the direction of the vector k is not fixed, this condition can 
always be satisfied in an infinite system for sufficiently longwave fluctuations. In polarized quantum gases, with the 
exception of small polarizations, the spin rotation parameter \iM is not small, /j,M > 1. Then the condition ( |l7j ) is 
roughly equivalent to 

k 2 < 2^ik • VM . (18) 

This type of instability in transverse channel was observed in several experiments with 3 He |, Refs. [|l2]-jl4j. 
Though the conditions and the setups varied, the experiments were performed for the systems in which the coherent 
effects exceeded the diffusion damping, /iMq > 1. Usually, the initial condition corresponded to the large gradient of 
longitudinal magnetization that could be created, for example, by a 180° pulse in a part of the cell. The instability 
started in accordance with Castaing criterion and often led to long-lived NMR ringing. Though the Castaing instability 
in transverse channel often serves as an initial stage of the domain formation, the instability, by itself, is not necessary 
for the domain formation (HI . 

The experiment with trapped 87 i?6 Q was performed in such a way so that the Castaing instability was not 
displayed. First, the polarization in the initial state was already transverse thus eliminating the need for strong 
transverse currents on the initial stage of the domain formation. Second, this transverse magnetization, probably, 
did not have large non-equilibrium gradients. However, the Castaing instability is within an easy reach and can be 
observed in a slightly modified experiment. 

In JILA-type experiments, the analog of the spin rotation parameter [iMq is TgM$. Numerically this parameter 
is /jMo = rgMo ~ 17A, where A = Mq/ti is the degree of pseudo-spin polarization. For trapped alkali gases, it is 
difficult to achieve a high degree of the pseudo-spin polarization using the so-called brute force technique, i.e., the 
effective "magnetic field" f2 l- On the other hand, the preparation of the gas with a high non-equilibrium polarization 
is even simpler than in helium and does not require elaborate techniques such as optical pumping, Castaing-Nozieres 
method, etc. For example, one can simply start from a pure component. 

The instability should be re-examined for trapped alkali gases for which the confining potential U imposes an 
additional limitation. The reason is that in a system of the size L, the wave vectors smaller than 1/L are absent 
and the condition ( |l7j ) must be considered in conjunction with the requirement k > 1/L. On the other hand, the 
polarization gradient VMo introduces another length scale, Lm ~ Mq/'V Mq, which, depending on the magnitude 
of the gradient and the length over which it is spread, can be larger or smaller than L. As a result, the instability 
condition in finite traps acquires the following threshold with respect to the trap size: 

"'^W (19) 

Roughly, the range of the " dangerous" wave vectors is 1/L < k < 2/zMo / Lm ■ 

For the instability to develop, the diffusion equilibration of the gradients should not be much faster than the time 
tc on which the instability develops. The diffusion relaxation of a longitudinal gradient VMo takes place on a time 
scale of the longitudinal spin diffusion. The lower limit for this time can be estimated as tp ~ L 2 /D. In certain 
situations of open systems the non-equilibrium gradients can be supported even on much longer times. Assuming 
that /iMq is large and k is far above the threshold, the estimate of the time tc ~ l/|Imw| from Eq.(|l6|) yields 

t c ~ L 2 M /D. (20) 

Comparison with the diffusion time shows that the instability in small traps can be observed only if 

L > L M . (21) 



4 



This is the simplest condition for the observation of the instability in finite isolated traps. 

If L ~ Lm, the diffusion time is comparable to the time of development of instability. However, since at /j,M ^> 1 
the instability is developed at relatively large wave vectors k > 1/L, the signs of the instability should be still clearly 
visible in experiment as higher harmonics in magnetization distribution in transverse channel, i.e., "ringing" similar 
to that observed in helium experiments, Refs. |T^-|l4|]. 

Above we neglected the possibility of having gradients of magnetic held which can stabilize the non-uniform dis- 
tribution of magnetization P,|fO|. In order to observe the instability in the presence of the held gradient, Eq.(pl|) for 
Lm should be supplemented by the condition that Lm is also much smaller than the width of the quasi-stationary 
domain wall A, Eq. (p^), 



L M < 



( D ) 1/3 ~L 


r 1 u:l 


1/3 


\fjLM\vn L \J 


AiranL 2 u 2 





(22) 



According to the estimates from the previous section, A > L/3 for the held gradients v < 0.1 Hz. This means that 
if the held gradients are noticeably smaller than this value, Eq. (p2] ) does not impose any new restrictions for the 
observation of the instability. 

Note, that in the situations, when the non-equilibrium gradient is supported externally, i.e., when the quasi-steady 
currents can how, the condition ( pll) is not required. Under such conditions, the scale Lm is Lm ~ L and the 
instability develops, as in helium systems, as long as fiMo 3> 1. 

There are many experimental ways for creating a non-equilibrium gradient of polarization of the trapped gases 
sufficient for the observation the Castaing instability, Eq. (^l|). As a result of this instability, the polarized system 
should acquire a rapidly growing perpendicular component, M fr ±M . The dynamics of similar processes in helium, 
though in very different geometries, has been studied computationally in Ref. ]To} ] . 

In experiment [0, the 90° r/-pulse created the almost uniform transverse polarization Mo in, i.e., the x-direction 
in the pseudo-spin space. The spin-rotation parameter was relatively large, 



liM » rgn/2 ~ 9. 



(23) 



[The trapping potential was responsible for some non-uniformity of polarization. This equilibrium non-uniformity 
could not lead to any instability; actually, small held gradients help to stabilize the system]. The small fluctuations 
(5M around this configuration were, if one ignores the small magnetic held, the gapless Silin waves, Eq. (|l^), and 
there was no instability. However, if the initial polarization M x in experiment (j were sufficiently non-uniform for 
the instability to occur, the observed component separation , i.e., the formation of domains with different M z would 
have taken place much faster than in Refs. Jl|-^|]. This can be checked experimentally. 

There are several other possible experiments with the instability and domain formation that can be performed with 
trapped alkali gases along the lines of the helium experiments. First, one can start from a one-component gas and 
flip the "magnetization" by 180°, i.e., excite the Rabi transition, in a part of a trap. Then one can observe long time 
magnetization ringing in the system similar to Refs. |l2|-]l4|]. In this experiment, however, the transverse polarization 
might already exist after the pulse on the interface between the components due to some residual non-uniformity of 
the Rabi frequency ]l5| ]. In this sense, it is not clear whether the instability can be unambiguously identihed. 

The next experiment involves bringing together two clouds consisting of the two different components, which have 
no memory of each other. Employing the pseudo-spin language, this two clouds are characterized by the opposite 
polarizations in the z-direction with no transverse components. Then, the dynamics of mixing of this two clouds can 
proceed along the line of Castaing instability, as opposed to a trivial (incoherent) diffusion mixing. Accordingly, the 
transverse non-uniform magnetization will arise spontaneously near the interface between the mixing components. 
Then it will slowly relax till the equilibrium state is reached. In some sense, this effect is the time reversal analog of 
the experiment [jjj . Indeed, in the case [jlj , the initial non-equilibrium transverse polarization resulted in the transient 
separation of the components. The effect of Castaing instability results in the spontaneous creation of the transverse 
polarization in the process of the components mixing. 

Another possibility is to create a non-uniform stable distribution of longitudinal polarization in which the criterion 
( pl| ) is violated and then to change rapidly the size of the trap so that to comply with the instability condition (|2l|). 

The last suggestion is more speculative. The non-uniformity of the initial non-equilibrium conhguration can be 
caused by the sensitivity of the two-photon Rabi transition to the gradients of difference of the trapping potentials 
p5[ ; that is, to the value of the Larmor frequency f^. This spacial non- uniformity of the pulse frequency f2 (r) 
results in formation of helicoidal structures after the pulse. For example, if the Rabi rotation takes place around the 
a;-axis, the structure of magnetization M (r) after the pulse of duration t has the form M z + iM y = Mq exp(i£t(r)t) . 
Such helicoidal polarization M (r) has large gradients and may be unstable. As was emphasized in Ref. |15|] , the 
Rabi frequency gradient Vf2 can be controlled with high accuracy. Accordingly, one can control the magnetization 
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gradients. However, the time evolution of the helicoid-like structure M (r) is described by a convoluted set of coupled 
transverse and longitudinal equations. At present, it is not clear to what extent the instability condition for such 
helicoids differs from the above equations for simple geometries. This general suggestion can be modified in various 
ways so that to create different structures some of which can be even one dimensional and are easier to describe 
theoretically than the helicoids of a general form. However, a practical implementation strongly depends on an 
available experimental setup. 

In most of the suggested experiments, the detection of the instability and the spontaneous formation of the transverse 
magnetization, i.e., of the large off-diagonal components of the density matrix associated with continuous Rabi 
transitions, can be done using the Ramsey spectroscopy Jl6[ |. 



V. CONCLUSIONS 



In conclusion, we pursued the analogy between the two-component Rb gases and spin-polarized helium systems. 
The observed separation of components in experiment Jl]] is analogous to the formation of coherent precessing domains 
in helium ||. This interpretation is consistent with the conclusions of Refs. [^]-[|]. The width of the domain wall 
agrees with the formation of up to three domains. An increase in the gradient of the pseudo-magnetic field can 
result in a spontaneous stratification of the system. A polarization gradient in the initial state can drastically speed 
up the domain formation. Apart from the observation of the domain formation, there are several other interesting 
experimental possibilities associated with the pseudo-spin Castaing instability in transverse channel. This instability 
is almost inevitable in highly polarized systems with non-uniform distribution of the pseudo-spin and can be observed 
in the JILA-type experimental setups. The conditions for this instability in finite traps are delineated. 

The authors acknowledge support by the CUNY grant PSC-63499-0032 (A.K.) and by the NSF grant DMR-0077266 
(A.M.). 
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